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imODUCTION 


Theoretical Inveetlgation of the breakdown of an laolated 
vortex haa bean conducted before (Rafa. 1 and 2). However, ita 
application to alender-wing aerodynaaica is quite limited in 
scope and versatility (Ref. 3). In fact, very few correlations 
with dsta have been reported. In addition, these available 
theoretical methods require lengthy computer time and hence, 
sre not suitable for applications in a preliminary design. 

In this paper, a saml-emplrical method to predict effects 
of vortex breakdown cn aerodynamic chr.racterlstlcs of slender 
wings is described. The method is based on Polhamus* method 
of suction anslogy. Both longitudinsl and lateral-directional 
characteristics can be predicted. 

THEORETICAL DEVELOPMENT 


In the method of suction analogy, the vortex lift is equated 
to the edge suction force in the attached flow. It was observed 
by Lamar (Ref. 4) that for delta wings the angles of attack for 
vortex breakdown at the trailing edge ( of geyg) can be related 
to the leading-edge suction disirlbution kC^C. ). For a large 
o<apTtf • " distribution tends to peak out at a more 

outboard location (Fig. 1). However, the location of ( Cf C 
itself is not a good correlation parameter because for lower swept 
wings, such as the 60-deg. delta, th peak location is less 
distinguished. 
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Tharcfort, th« method to be developed muet (I) employ e good 
correletion peremeter for eppllceble to moet plenforme, 

not Juet delta vinge; (2) be able to predict the progreeslon of 
the breakdowt point when •tytff.gprt • predict the amount of 

remaining vortex lift in the region of breakdown. Theee are 
described in the following for both longitudinal and lateral- 
directional aerodynamics. 



where c i> the mean aerodynamic chord. Since Cg Is proportional to 
5/eV . it follows that ^ is a function of planform and Mach 
number only. If ft is the nondimenslonal centroid location of the 
^ - distribution from inboard to f ”** found that 

^ correlated well with experimental oipptt shown In Fig. 2. 
The relation between and Iff la illustrated in Fig. 1. The 
experimental are taken from Reference 5. The correlation 

curve in Fig. 2 was obtained by a least-square analysis of Wentz 
data (Ref. 5) and can be described as follows. 


^»»TB 


tnt - 




(2b) 


In Reference 5, the progressirn rate of the breakdosm point 
on delta wings was also measured. Note that the progression rate 



Masurcd in th« water tunnal la aoMwhat diffarant from that in 
tha wind tunnal (Raf. 6). Although tha axparlaaintal prograaaion 
rata la alightly diffarant for diffarant delta wlnga (Fig. 3) , for 
alnpliclty, a alngle curve obtained by a laaat*aquara analyaia of 
Wanta data ia uaad in tha praaant analyaia. Tha curve can ba 
deacribad by 

^ ;/ «a( < AO Wy. (3a) 

AX 3 ^ (3b) 

where a«(« (4) 

and AX la the nondinenaional It “ diatance (referred to the root 
cho'id) from the trailing edge to the breakdo«m point. The apanwiae 
location of the breakdown point, , nondlmenaionalized with half 

spsn, can therefore be obtained as 

Once the vortex - breakdown point ia located on the planform, 
it la necessary to determine the amount of remaining vortex lift 
after breakdown. It is wadi known that the vortex strength after 
breakdown is decreased, but not vanished. In other words, the 
sectional at any spanwlse station y ^ anist be 

t&ultlplled by a factor 4 which is less than 1.0. To determine 4 , 
Wentz data for delta wings were analyzed at A«<m^deg. and compared 
with results by the Quasi-Vortex-Lattlce Method (QVLM) (Ref. 7). 

The factor 4 is determined by requiring the QVLM results to match 
Wentz data as closely as possible through a least-square analysis. 
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Tha raaults ara 

4 m •./!/ -h J 

4 a -tf#/- ?/ (6) 

Tha raaulta indicata that dlffarant planforms (i.a. diffarant fjf ) 
hava dlffarant vortax charactarlatlca aftar braakdown. 

Lataral-Dlrectlonal Aarodynaalca 

It waa shown In Rafa. 8 and 9 that In a sideslip, tha laading- 
adge vortex on the windward side Is pushed Inward and downward, 
thus Inducing additional vortex lift; while on the leeward side. It 
Is pushed outward and upward, and thus Inducing less vortex lift. 

Note that the usual vortex lift, as explained through the suction 
analogy. Is produced by the existence of net upwash at the leading 
edge, thereby Inducing vortex separation. This Is called the roll-up 
type of vortex lift in the present analysis. On the other hand, 
at a given ot the sideslip does not change the net upwash at the 
leading edge. Therefore, to account for the Increase In vortex lift 
on Che windward side due to sideslip, a "dlsplacenenc-type" vortex 
lift must be Introduced. The flow tnechanisns and computational 
methods for sideslip, yawing and rolling motions are described below. 

(1) Vortex lift of the displacement type in sideslip 

Tn general, the leading-edge thrust coefficient is given by (Ref. 10) 

g. . j. (7) 

where ^ Is Che leading-edge sweep angle, and C is the leading-edge 
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singularity paramatar daflnad as 


c 


Litn 




( 8 ) 


Now, In sldasllp, Is glvan by (Raf. 11) 

whars the t -coordinate In the fraaatraan direction and ^ 

are the components of vortex density In the X - and / -> directions , 
respectively. The component, , Is responsible for the longitudinal 
loading. Since the vortlclty vector near the leading edge must be 
parallel to the leading edge. It follows that 


^ 3|f Aj 


( 10 ) 


near the leading edge. Hence, by substituting Eqs. (9) and (10) Into 
(8) , It is obtained that 


C- CU(/± UnAiHof) 


(U) 


where 



( 12 ) 


is the leading edge singularity parameter due to of for the symmetrical 
loading. Eq. (11) shows that the singularity parameter due to sideslip 
can also be Interpreted as that produced by the Interaction between 
and the sideslip velocity at the leading edge: 



Lm 



(13) 
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By squarlog Eq. (11) and aubatitutint it into Bq. (7), tha laadlnt- 
adga thruat coafflciant in aldaallp can ba obtalnadt 


net + *«>% y*V)V/-MW4< 

(14) 

Tha raaulting thruat la In tYm dlractlon pointing Into tha aldaallp 
valoclty (Fig. 4a). Ranca, tha corraapondlng auction ccafflclant 

CL « - Ct f t 0 

^tf) fi»Af \ r ' X f/ 


who re 




. 


/7 




(16) 


la the auction coefficient for the aynmetrical loading. 

(2) Vortex lift of the dlaplacement type In yawing 

Aa ahown in Fig. 4b, the Induced aidewaah at la 



Vi 







(17) 


Thla aidewaah varlea over the planfom, and hence givea rlae to 
"variable aidealip" effect. According to Eq. (15), the correaponding 
auction coefficient ( ) la then 




(18) 


(3) Regiona without potential flow effect 

In the method of auction analogy, the vortex lift la added to 
the potential-flow lift to obtain the total lift. The potential- 
flow lift ia calculated by aumning all attached-flow preaaure forcea 
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ov«T the planform. This is uasd for Che situation in which a vortax 
flow of cha roll-up type is involved, such as in tae syoBstrical 
loading or due to rolling aotion. For a vortex flow of the dis- 
placeaent type, a certain region near the edges has to be excluded 
in calculating the potential-flow cosiponent to avoid over-prediction. 

According to Reference 12, the magnitude of gives a good 

estimate of the location of vortex lift centroid from the leading 
edge. Hence, ^ and are taken to be the streamwlse 

distances measured from the leading edge to the vortex lift centroids 
for ^ and Y nodes of motion, respectively. In the present method, 
regions near the leading edge with ^ are assumed 

not to nave the potential-flow effect. However, if the vortex 
breakdown is predicted to occur, the vortex flow region tends to 
enlarge, so that the symmetrical ( ^ C ) will be used. 

(4) Boundary layer effect on Roll Damping 

The rolling motion induces a change in local angle of attack 
so that the roll-up type vortex lift will be produced. Its calculation 
is described in Reference 7. 

However, It Is known that the conventional attached-flow 
theory tends to over-predlct the roll damping by as much as 20X 
without accounting for the vortex flow (Ref. 13). One possible 
reason Is that the boundary layer near the wing tip is thickened 
by the centrifugal effect in rolling, so that the effective angle 
of attack is reduced. For engineering applications, this effect 
will be accounted for in the following way. 

Note that the boundary layer thickness is a function of skin 
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friction co«£ficl«nt ( ). In R«f«r«nc* 14, a foroula for tha 

airfoil lift raduction dua to boundary layar thlcknaaa la givan aa 




iV9 


(19) 


whare 4' choaan as 2.0. For a wing in staady rolling, it la 
asaused that Eq. (19) la appllcabla at tha station of naan aarodynaolc 
chord (MAC). According to Rafarenca 15 (page 599), tha thickness of 
a 2-D turbulent boundary layer can be wrlttei. as 

( 20 ) 


where Tff indicates that is to be calculated at the trailing 
edge and Jff is the kinematic viscosity. It follows that the slope 
of boundary layer is given by 


Jf, 




__ y ly 

0.37 ^ 


( 21 ) 


According to Ref. 15(page 607), the boundary layer thickness on a 
rotating disk is given by 


* o. itu 


where ^ in the present application is taken to be the y - coordinate 
of MAC and cO is the rotating speed. It follows that 





(23) 


As sho%m in Fig. 5, in the present 3-D applications, 
is taken normal to the trailing edge. The effect of steady rolling 
is regarded as a small perturbation to the angle-of-attack effect, 
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so that only its component along the trailing edge will be accounted 
for. The angle, , for the direction of combined slopts Is 
then given by 



X -f~ 



(24) 


However, 


CC 





(25) 


Hence , 

V^Xt V^Xt Xt L 

Substituting Eq. (26) into Eq. (24), it is obtained that 

tia>» = f.»46X. cos 


(26) 


(27) 


where 


(28) 

It Is assumed that if Eq. (19) is fcr the angle-of-attack 
effect, then the combined effect , given by Eq. (19), is increased 
by a factor of l/C 0 Sf^ with given by Eq. (27). With the local 
angle of attack due to steady rolling, , replacing ^ in Eq. (19). 
Eq. (19) is now rewritten as 

* 7 / (29) 

where the constant is taken as 2.2, instead of 2.0, for a better 
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ccrrclation wlch data, and 1* Included for tha notching 

affect of the trailing edge. That ie, the component, » 

shown in Fig. 5 will make th»'> traillng-edge boundary layer thinner. 
(5) Rolling and yawing at non-zero sic slip 

It was shown in Reference 8 that the roll damping of a slender 
wing at high angles of attack car be Increased if the wing rolls 
about the stability axes at a non-zero sideslip. To explain this 
effect, let the total leading-edge singularity parameter be 

C * C^(! ± 'Un f ) + ^ ^ Ai) **■ ^ 

where and are leading-edge sing;.0 arity parameters 

produced by upwash (i.e. the roll-up effect), respecf’vely. By 
squaring C to be used in Eq. (7), various cross-product terms appear. 
For the rolling motion, the vortex lift due to roll-up is further 
Increased on the right wing due to Inboard displacement of vortex 
core in positive sideslip. Similarly, any roll-up effect in the 
yawing motion \t 111 be further Increased or decreased due to sideslip. 
However, products of any two roll-up effects or displacement effects 
do not produce forces and moments of the antlsysnetrical nature. 
Hence, additional vortex-lift terms in the leading-edge singularity 
parameter due to rolling and yawing in sideslip are 

^ ^^Aj f) (31) 

Eqs. (31) and (32) are to rep ice in Eq. (7) to produce the 
additional leading-edge thrust. 
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(6) Vortex br*akdo%m in sideslip 

In general, the angle of attack for vortex breakdown at the 
trailing edge in sideslip, , is reduced on the windward 

side and increased on the leoward side. The following steps are 
used to calculate present analysis. 

(a) Calculate ^ froa the syaaMtr''cal suction distribution 

( Cg ) , except that the leading-edge sweep angle 

is taken to ^e Aj-f for Ch. right wing ^ 

for the left wing. 

(b) Assuae that the aaxlaua vortex strength before breakdown 
for a given planfora is unchanged by sideslip. Since 
the vortex strength is represented by ^ which is 
proportional to sln^e< , it follows that the aaxlaua 
vortex strength will be reached at a lower a< ' than 0(g^jg 
on the right wing because of Increase in vortex lift 

in sideslip. Hence, can be sclved from the 

following equation 

(c) The vortex lift factor ( ) is calculated by the saae 

formula for (i.e. Eq. (6)), except that is 

replaced by calculated in step (a). 

(d) Whenever ihe vortex breakdown is predicted to occur at 

a given y station. itt) <»ust be multiplied by 

Therefore, the antisyooetrical vortex lift on the windward 
side due to sideslip would be given by 
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On the leeward side, the vortex breakdown may not occur at a given . 

In this case, the magnitude of the antia 3 nnmetrical vortex lift is the 
average of those on both sides. 

(7) Vortex breakdown in yawing 

The effect of yawing on vortex breakdown is similar to that due 
to sideslip. With a positive yawing, the left tip will be the first 
which is subject to possible vortex breakdown effect. Since the equi- 
valent sideslip is variable along the leading edge, the average of 

at a given y station and that at the tip is used to determine 
For simplicity, for the yeving motion is determined by a linear 

interpolation between ^BDTE sideslip . 

When the yawing motion occurs at f 0, the equivalent sideslip to be 
used should be -f-fj. on the left wing and^ + ^^ on the right wing. Since ^ 
is usually negative near the tips, the equivalent sideslip there is con- 
siderably reduced. In this case, the vortex breakdown may not occur, 
in particular at a high yaw rate, so that the rolling moment due to 
yaw rate may be large. 

(8) Vortex breakdown in rolling 

Due to a positive roll, the local angle of attack on the right 
wing is increased. It would seem that the vortex-breakdown of would 
be decreased. However, rolling tends to move outboard the centroid of 
the vortex lift distribution (i.e. increasing y^) , so that the vortex- 
breakdown of is increased. In the present theory, it is assumed that 
for the rolling motion is equal to for the symmetrical loading. 

The average of the local of at a given y station and that at the 
tip is used to determine whether there is vortex breakdown at that 
station. 
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(9) Cancroids of vortex lift 


Based on the Idea developed in Reference 12, c c gives the location 

8 

of vortex lift centroid. Let leading-edge suction coefficients produced 
by upwash at the leading edge be denoted by *^s(r)' 

roll-up for the antlsytnBetrIcal vortex lift is assumed to cake place along 
a direction perpendicular to Che leading edge. Hence, Che antlsymmecrlcal 
vortex lift of the roll-up type is assumed to act at 

/ - I 1 Ajt 

for the rolling motion, and at 


for the yawing motion. In sideslip, it is at 
This is illustrated in Fig. 6a. 

On Che other hand, if rolling and yawing motions take place in a 
non-zero sideslip, the vortex lift centroids will be at 


y- 



SfH Ax 


y- 






Aj 


For rolling and vawing motions, the Inboard displacement of vortex 

lift centroids is further restricted by the position of symmetrical vortex 

lift centroids (i.e. R ■ c c) . Therefore, if 

s 

l<^,c I f/»Aj > At As 

it is replaced by R cos /l^sinyl^, where R cannot exceed the local chord. 
This is illustrated in Fig. 6b. 

In a coning motion, rolling and yawing motions are coupled. The 
vortex lift centroid Is then assumed to be that due to the rolling 


mot ion. 
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NUMERICAL RESULTS 


Longitudinal Characteristics 

Results for delta wings of 60*^ , 70° and 80° are presented in Figs. 7- 
10. The method for computing the vortex lift involving rounded leading 
edges is based on Kulfan's concept (Ref. 16). Note the differences in 
lift curves after vortex breakdown for different delta wings. The 
lift curve for the 80-deg. delta has a sharper change after vortex 
breakdown as compared with that for the 60-deg. delta. The lift curve 
for the 70-deg. delta wing represents the transition between these two 


types. 

Lateral-Directional Characteristics 

The effect of sideslip on the angle of attack for vortex breakdown 
at the trailing edge is presented in Fig. 11. The trend is correctly 
predicted, although the predicted magnitude is 2 to 3 deg. off. 

The lateral characteristics for 60-deg. and 70-deg. delta wings 
are presented in Figs. 12 and 13, respectively. The predicted results 
agree well with data, except when of exceeds by 8 to lU deg. In 

the latter situation, the effect of boundary layer separation must be 
accounted for and is not included in the present method. 

Additional results are compared with data in Figs. 14, 15 and 16. 
As shown in Fig. 15, the boundsT'y layer correction is seen to be impor- 
tant for accurate prediction of roll damping. The change in slopes of 


C. and C- after vortex breakdown is due to the difference in breakdown 
•t's for left and right wings. 

Finally, the effect of sideslip on roll damping is illustrated in 

Fig. 17. With zero sideslip, the slope of -p curve is always positive. 

This means that C is positive. On the other hand, as the sideslip is 
'*P 

increased in a steady roll, C. may become negative, as shown by the curve 




representing / -10°. At|^-5° and small p, there is some discrepancy 

between the theory and data. 

Extension of the theory to other planform shapes is currently being 
investigated. 
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